Role of anisotropy in the Forster energy transfer from a semiconductor quantum well 

to an organic crystalline overlayer 
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We consider the non-radiative resonant energy transfer from a two-dimensional Wannier exciton 
(donor) to a Frenkel exciton of a molecular crystal overlayer (acceptor). We characterize the effect 
of the optical anisotropy of the organic subsystem on this process. Using realistic values of material 
parameters, we show that it is possible to change the transfer rate within typically a factor of two 
depending on the orientation of the crystalline overlayer. The resonant matching of donor and 
acceptor energies is also partly tunable via the organic crystal orientation. 
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I. INTRODUCTION 

A large effort has been devoted to the successful devel- 
opment of organic light emitting diodes (OLED) despite 
their intrinsic limitations, particularly regarding carrier 
injection and transport, as compared to inorganic semi- 
conductor devices. Inorganic-organic hybrid materials 
may allovi' the rules that apply to matter and light to 
be further stretched, potentially seeding a new paradigm 
in optoelectronic devices taking advantage of the best 
of both worlds^. This is of paramount relevance in the 
strong coupling regime in which novel hybrid quasiparti- 
cles are formed, but also, in the weak coupling regime in 
which Wannier and Frenkel excitons maintain their indi- 
viduality, an hybrid system offers significant advantages^. 
In particular, a way to circumvent the drawbacks of or- 
ganic materials mentioned above is to use an heterostruc- 
ture containing an inorganic semiconductor subsystem in 
which carriers are electrically injected, transported and 
bound into excitons, coupled to an organic light emit- 
ting subsystem via a Forster energy transfer proces^. 
As a step in this direction, it has been proposed and re- 
cently demonstrated that non-radiative energy transfer 
can be efficient enough (see Refi for a review of relevant 
work) . Following early theoretical predictions^"—, the en- 
ergy transfer process has been observed from a quantum 
well to a quantum dot overlayeri^^, and from a quantum 
well to an organic overlayer, see Illl - KLSl to name a few. 

We deal here with such hybrid systems structured in a 
planar geometry whereby a Wannier exciton in an in- 
organic semiconductor quantum well plays the role of 
the donor and a Frenkel exciton in a crystalline or- 
ganic overlayer that of the acceptor. While many or- 
ganic subsystems used so far are not crystalline and ef- 
fectively isotropic, most of the organic materials of in- 
terest are strongly anisotropic if grown as oriented single 
crystals. Previous theoretical work has only considered 
the isotropic case and here we extend it to include the 
optical anisotropy of the organic subsystem. This will 
provide a complete quantitative description for the case 
in which the organic acceptor layer is an oriented sin- 
gle crystal and will allow to characterize how the energy 



transfer process depends on the orientation of its princi- 
pal axes. In section |lll we describe the theory of Forster 
energy transfer in a planar hybrid nanostructure taking 
into due account optical anisotropy. In section IIIIl we 
employ a model optical dielectric tensor and parameter 
values of typical organic media to calculate the transfer 
rate for various configurations, and we discuss the role of 
anisotropy. 



II. THEORETICAL MODEL 

We consider here the planar architecture shown in 
FiglH The donor subsystem consists of a semiconduc- 
tor quantum well of thickness 21 sandwiched between two 
semiconductor barriers of thickness {V — I). For simplic- 
ity, we take the same background dielectric constant Sb 
for the well and the barriers and we assume the barriers 
to be infinitely high (i.e., the Wannier exciton is fully 
contained in the well region). At the bottom [z < —I') 
lies a transparent glass substrate with dielectric constant 
Eg while at the top (z > I') lies the acceptor subsystem 
consisting of a crystalline organic medium with dielec- 
tric constant Sij . Both of them are supposed to be semi- 
infinite. The quantities Sb and Eg include only the contri- 
bution of higher resonances (with respect to the exciton 
energies under consideration) and we consider them to be 
real. The quantity is the total dielectric function of 
the anisotropic organic material, including in particular 
the resonant absorption due to the Frenkel excitons, and 
thus it is a complex valued tensor. 

Here we use the same theoretical framework to calcu- 
late the Forster energy transfer rate in a hybrid nanos- 
tructure as in references Hi, and discussed in detail in 
reference 0. This approach is equivalent to the usual 
Forster theory based on the dipole-dipole interaction— 
and leads to a macroscopic, semiclassical description of 
the energy transfer. In short, the transfer rate is obtained 
from the Joule losses suffered in the organic medium by 
the electric field induced by the quantum well exciton. 
As a matter of fact, the presence of the Wannier exciton 
of energy Ee^c = ^ gives rise to a source term in the 
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j{z) = EbSij at —I' < z < V 



and anisotropic (and complex) in the organic layer: at 
z > r, Sij{z) — Sij. The appropriate boundary condi- 
tions at the interfaces are the continuity of the tangen- 
tial electric field and of the normal electric displacement 
between each layer. Knowing the electric field we can 
calculate the transfer rate (inverse transfer time) of the 
excitation into the organic medium: 



FIG. 1: Sketch of the planar hybrid heterostructure. The 
inset shows the angles used to define the orientation of the 
organic crystal in the top layer (see text for details). (Color 
online) 
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quantum well corresponding to the exciton polarization 
oscillating at frequency lo which can be written as: 

P(r)e-*"* = d''=V(i-,i-)e"'"*, (1) 

where d"'^ is the matrix element of the electric dipole 
moment between the Bloch functions of the conduction 
and valence band extrema and ';/'(re,r/i) is the envelope 
function describing the bound electron- hole pair (re, Yh 
being the electron and hole coordinates). In our simple 
model of a 2D Wannier-Mott exciton, the polarization is 
given by the product of the Is-wave function of the 2D 
relative motion of the electron-hole pair with the lowest 
subband envelope functions of electron and hole (which 
are equal in the approximation of infinitely deep well) 
and with the plane-wave function of the center-of-mass 
motion. Taking into account their corresponding normal- 
ization we can write: 

where S is the in-plane normalization area, K the in- 
plane wave vector of the center-of-mass motion, r|| = 
(x,?/) the in-plane component of r and is the 2D ex- 
citon Bohr radius. 

Then we look for the electric field E(r) g^*"^* resulting 
from this polarization. Having in mind a distribution of 
Wannier excitons having a broad range of wavevectors as 
obtained by non resonant pumping or electrical injection, 
or simply a thermal distribution corresponding to a tem- 
perature of the order of 100 K, we can neglect retardation 
effects in Maxwell's equations {k >> uj/c limit) and use 
an electrostatic approximation^. The electric field can be 
written as E(r) = — V(/3(r) and is found from the Poisson 
equation for the potential 

V,ey(z)Vj(^(r) = -47rp(r) (3) 

where p(r) = —V • P(r) is the source charge density, and 
all the monochromatic oscillating factors e~'"* will not be 
explicitly indicated. The dielectric tensor £ij{z) is piece- 
wise constant corresponding to each different layer. We 
choose it to be isotropic (and real) in the inorganic part: 



This expression is equivalent to applying the Fermi 
Golden Rule to the decay of one excited state in the quan- 
tum well into the excited states of the organic molecules 
in the linear regime approximation (for a derivation of 
EqlHin terms of a microscopic inelastic scattering rate of 
Wannier excitons due to resonant two-level molecules in 
the organic medium see Reffisl). The power W dissipated 
in the organic layer is given hy W = Huj/t. Such energy 
transfer mechanism have been shown to be fast enough to 
efficiently quench the Wannier exciton luminescence and 
to turn on the organic molecule light emission. While pre- 
vious theoretical work has only considered the case of an 
effectively isotropic organic medium (i.e., iij ex. Sij), we 
focus here on the effects of the anisotropy of the organic 
layer on the transfer time r, assuming for the donor sub- 
system the same model appropriate to zincblende semi- 
conductors used therei. 

Due to the in-plane translational symmetry of the 
source (we are dealing with free excitons in the well), 
we consider the polarization for a given in-plane wave 
vector. In this case three modes of different symmetry 
can be identified: longitudinal (L) where d"^ is along the 
in-plane wave vector, perpendicular (Z) where the dipole 
moment is oriented along the z-axis, and transverse (T) 
for which the polarization is orthogonal to the two first 
case. We notice that the T mode does not give rise to 
any charge density and thus the dipole-dipole interaction 
vanishes. The two remaining polarizations lead respec- 
tively to two different charge densities in the quantum 
well p{r) = p(^'^)(z)e*^'"ii , given by: 

p^^^z) = -iKlpoil + cosqz) (5) 
p^^'' (z) ^ qlpo sin qz (6) 



For any given source exciton symmetry (L or Z, which 
will not be explicitly indicated), writing ip{r) = (/)(z)e*^'"ii 
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the Poisson equation gives 



plane wave vector 'K = K (cos A, sin A) is given by 



dz^ 
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In the organic part (f) has the form — poC e'*'^*^^ ' •* 
with C an amplitude coefficient and 7 a complex number. 
The value of 7 is the solution of the last equation in 
such that i?e(7) < to satisfy the boundary condition at 
infinity: 
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consequently 7 depends on K and on the orientation of 
the organic material. Solving the 8 by 8 system obtained 
imposing the two boundary conditions on each of the 
four interfaces (at z — ±1 and z — ±{l + 1')), we find for 
the two source modes, respectively, the amplitude coeffi- 
cients: 



(L) ^ 

sinh(XO {ei, cosh{Kl') + Eg sinhiKl')) 
{Sg + e)£b cosh(2Kl') + {sgi + el) amh{2Kl') 
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where the parameter e 



-£zzl 



• K 



- Ky 



comes from the continuity of the normal electric displace- 
ment at the boundary between the organic and inorganic 
layers. These amplitude coefficients depend on K and 
on the orientation of the organic material in a rather in- 
volved way through e and 7. We recover previous results 
for the isotropic case^ provided that lij oc Sij and eg — e. 
The electric field in the organic layer is then: 



E(r) = HK-7^e,](^(r) 



t{K,X) 2n^ha% {21)^ Rej 

Imexx cos^ A + Ime^y sin^ A + ImCxx I7P 
+ 2Im7(Ime2,2 cos A -I- Ims^z sin A) 



2Im£a;^ cos A sin A 



(13) 



(12) 



The increased complexity of the present analytical re- 
sults with respect to previous theoretical work! is due 
partly to having lifted the mirror symmetry with respect 
to the z = plane assumed there, and partly to having 
included an anisotropic dielectric tensor for the accep- 
tor subsystem. The latter point gives rise to significant 
consequences on the energy transfer process as discussed 
below, while the former point allows to describe the quan- 
titative dependence of the transfer time on the value of 
the substrate dielectric constant which, however, within a 
reasonable range of parameter values is a minor effect. In 
most cases of experimental interest, all directions for the 
center of mass wave vector of the exciton in the well are 
equiprobable (i.e., the exciton distribution has cylindri- 
cal symmetry), and an average over the angle A is taken, 
which is done numerically. For a thermalized population 
of quantum well excitons, we also average over the en- 
ergy , with M the exciton mass, according to the 
Boltzmann distribution. Taking these averages tends to 
reduce the effect of the anisotropy, and this would also be 
the case for the energy transfer of a localized exciton, in 
which case the wave vector distribution is given by the 
Fourier transform of the localised wave function of the 
center of mass motion. 



III. RESULTS AND DISCUSSION 

For illustrative purposes, we will focus on the effects of 
the anisotropy and thus only change the organic overlayer 
configuration keeping all the rest of the heterostructure 
fixed, and similar to that assumed in previous theoret- 
ical work. We will thus consider a variety of cases of 
anisotropic organic materials without bothering to select 
in each case an appropriately matched donor subsystem. 
In this way, the effects of the anisotropy will be singled 
out and discussed. For the donor subsystem, we take a 
II- VI semiconductor quantum well (e.g. ZnSe) for which 
typically-Si ~ 6, d^'^ ~ 12eaB and as is taken to be 
25A, M = 0.76too. For the structure, we take I — 30A 
and V = 40^, while the glass substrate and the organic 
top layer are infinite in the z-direction. 

The dielectric function of an anisotropic organic mate- 
rial can be written as^^: 



where is the unit vector in the z direction. Finally, 
from Eq. ([?]), the transfer time t for an exciton of in- 
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where £00 is the isotropic high-frequency dielectric con- 
stant, V is the volume of the lattice cell, d„_i is the zth 
component of the transition dipole moment of the nth 
excitonic eigenstate with energy fkUn and damping r„. 
In simple cases, it is possible to use an effective model^^, 
typically near the exciton resonances, where the optical 
response is modeled by a real background constant and 
several Lorentz transitions for each principal axis j, de- 
scribed by their energy Ej^, coupling amplitude Aj and 
damping Tj : 



iTjhu! 



(15) 



Ej^oo is here the background constant adapted to each 
axis. The values of the parameters are obtained by fit- 
ting experimental data. Here we use such a simple model 
with experimental values for the energy, coupling ampli- 
tude and damping of each relevant exciton transition, see 
Tables Hini 



TABLE I: Experimental values obtained for tetracene from 
the fitting of ellipsometry dataiS, valid for an energy range 
from 2.2eV to 2.5eV. 





El 


£2 


£3 


£00 


1.39 


1.00 


2.10 


Eq 


2.38 eV 


2.46 eV 




A 


2.971 


0.391 




r 


0.088 eV 


0.057 eV 





around followed by a rotation of angle ^' around u^: 

7^^7^"-(^-)7^^(0) 

This rotation can be rewritten in the device frame as 
n = W{e)n%'i'). Then the dielectric tensor e of the 
organic compound reads: 



(16) 



where £/ is given by (1151) . This orientation is chosen (or 
fixed) by the growth conditions of the organic compo- 
nent. In the following, we discuss the dependence of the 
transfer time on the orientation of the organic crystal, 
and the optimal choice of frequency for the donor for a 
given orientation. 




FIG. 2: Transfer time of a L-polarized exciton in function of 
the orientation for a uniaxial crystal with relative permittivity 
given by ei = £2 = 2, £3 = 2 -I- Si, dashed curve, and £1 = 
£2 = 2 -|- 3i, £3 — 2, solid curve. The straight lines are the 
mean values associated to the corresponding isotropic case. 



TABLE II: Experimental values obtained for a— tetraphenyl- 
butadiene (TPB) from the fitting of ellipsometry data^ ~ (con- 
verting the Gaussian fit into the Lorentz model, keeping the 
same peak value and area under the peak of the exciton res- 
onance), valid for an energy range from 3.2eV to 4.2eV. 





£1 


£2 


£3 


£00 


3.1 


2.84 


2.37 


Eq 


3.83 eV 


3.61 eV 




A 


2.46 


1.73 




r 


0.59 eV 


0.58 eV 





The orientation of the organic crystal is given by three 
Euler's angles (i", 9, \1/), as defined in Ref. [2(J : a rotation 
of angle $ around the z-axis followed by a rotation of 
angle 9 around the new x-axis and a rotation of angle ^' 
around the last new obtained z-axis. However, due to the 
rotational invariance of the inorganic subsystem around 
the z-axis, only the angles {9, 5') are needed. Besides, due 
to the mirror symmetry perpendicularly to each axis, it is 
enough to take them between and 7r/2. The coordinates 
of the optical axes {u^, i = x,y, z} expressed in the device 
frame {e,} (see FigH]) are given by = TZei where TZ 
is a rotation matrix constituted by a rotation of angle 9 



We consider a thermal distribution of excitons in the 
well and we assume a temperature of lOOK. At a given 
frequency of emission, we calculate the transfer time as 
a function of the orientation of the crystal. At first, we 
discuss the case of a uniaxial medium. Its orientation 
is determined by only one angle 9 as defined in p6p be- 
tween the optical axis and the normal of the planar struc- 
ture. In Figl2]we give the transfer time as a function of 
9 for a medium where the optical axis is the absorbing 
one (dashed curve) or the only non-absorbing one (solid 
curve) . We find a transfer time of the order of a few tens 
of picoseconds comparable to that of the isotropic case^. 
This is competitive to any other decay channel for the 
excitons in the well proving the efficiency of the energy 
transfer to the organic media. We observe, in particular, 
a variation of a factor 1.7 in the transfer time depending 
on 9, so the efficiency of the transfer is affected by the 
orientation of the crystal in a significant way. The mini- 
mum transfer time is lower in the case of two absorbing 
axis as expected as dissipation in the organic is increased 
in this case. 

To illustrate the case of biaxial media, we con- 
sider tetracene and Q;-l,4,4-tetraphenyl-l, 3-butadiene, 
for which dielectric function data were given above, and 
anthracene and PTCDA. For both the latter, we can 
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FIG. 3: Transfer time of a L-polarized exciton as a function 
of the orientation of the tetracene crystal, given by the two 
angles ^ and 6, at an energy 2.38 eV. for which e\ — 2.59 + 
3.0i, £2 = 2.21 + 0.05i, £3 = 2.09. (Color online) 




FIG. 6: Transfer time of a L-polarized exciton as a function of 
the orientation of anthracene crystal, given by the two angles 
and 6, at an energy 3.12 eV where ei = 4.71 + 4.86i, £2 ~ 
4.65 + 0.71i, £3 = 3.01 + 0.46i (ref. [U). (Color online) 




FIG. 4: Transfer time of a L-polarized exciton as a function of 
the orientation of the a-l,4,4-tetraphenyl-l, 3-butadiene crys- 
tal, given by the two angles ^' and 0, at an energy 3.83 eV 
for which £1 = 3.1 + 2.46i, £2 = 2.0 + 0.95i, £3 = 2.37. (Color 
online) 
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FIG. 5: Transfer time of a L-polarized exciton as a function 
of the orientation of PTCDA crystal, given by the two angles 
5' and 9, at an energy 2.2 eV. where £1 = 2.46 + 0.15z, £2 ~ 
7.1 + 3.41i, £3 = 6.3 + 6.5i (ref.iil). (Color online) 



use data obtained directly from ellipsometry mesurement 



at the relevant exciton energy in Ref. [2l|. We show in 
Fig. [21 m [S] and H] their respective transfer times as a 
function of the orientation of the crystal, defined by the 
two angles and 6 (ITBl) . The order ofmagnitude of the 
exciton lifetime is of a few tens of picoseconds, with a 
minimum around 20 ps. We observe a difference in the 
transfer time depending on the orientation of a factor 1.8 
for tetracene, 1.6 for TPB and PTCDA, and a factor 1.5 
for anthracene. The qualitative shape of these plots de- 
pends on the number of dissipative axis. In the case 
of a-1, 4, 4-tetraphenyl-l, 3-butadiene, the transfer time 
is higher when the non-absorbing axis is along the z- 
direction. For tetracene almost all the dissipation is along 
a given axis and the transfer time is maximum when one 
of the two non dissipative axis is along the z-direction. 
We show the one exciton life-time for a L-mode exciton. 
The Z-mode gives nearly the same results. 
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FIG. 7; Transfer time of a L-polarized exciton as a function 
of the resonant energy for q-1,4, 4-tetraphenyl-l, 3-butadiene 
crystal in two different orientation {6 = 7r/2,^' — n/2} thick 
curve, and {9 = n/2, ^ = 0} dashed curve. The minimums 
correspond to the two frequencies of absorption (cf Table |Il]) 
so it is possible to excite predominantly one or the other of the 
two excitons of the a- 1,4, 4-tetraphenyl-l, 3-butadiene crystal. 

These results have to be compared to the case where 
the organic compound is on a powder form, where the size 
of a mono-crystal is much smaller than the wavelength of 
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the resonant light. In this case an average over all possi- 
ble directions for the mono-crystals is realized, which is 
equivalent of averaging the dielectric tensor. We found a 
transfer time of 28ps for tetracene, 26ps for TPB, 21ps 
for anthracene and 17ps for PTCDA. Thus, in general, 
the transfer efficiency for this effectively isotropic case is 
comparable to that of a nearly optimally oriented single 
crystal. 




TV 



FIG. 8: Transfer time of a L-polarized exciton in function 
of the direction of this latter and of the orientation S of a 
uniaxial crystal, with relative permittivity given by ei = e2 ~ 
2, es = 2 + 3i. The transfer time is going to infinity for 
{A = 0, tt; e = 7r/2}. (Color online) 

For a given orientation of the crystal, as it could be 
imposed by the growth conditions, one can find a fre- 
quency that minimize the exciton life-time. In FigIT] 
we show the transfer time as a function of the reso- 



nant energy for two different orientation of the a-1,4,4- 
tetraphenyl-1, 3-butadiene crystal, each curve presenting 
a minimum respectively around 3.7eV and 3.9eV which 
correspond to the two frequencies of absorption of the a- 
1, 4, 4-tetraphenyl-l, 3-butadiene, each exciton being ex- 
cited separately. As a consequence, as a fucntion of the 
orientation, it possible to improve the transfer efficiency 
by improving the tuning with the emission frequency of 
the donor. 

Finally, for illustrative purposes, we consider the en- 
ergy transfer time given by Eqll3l for an excitons having 
a center of mass K along a preferential direction with 
a thermalised energy. We observe the variation of the 
transfer time with that direction given by the angle A 
and the orientation of the crystal 6. In the case of a uni- 
axial crystal with absorption in one direction, see Fig|8j 
there is no absorption when the Frenkel exciton is per- 
pendicular to the plane defined by the source exciton and 
the z-axis, then the transfer time is infinite. In the case 
of tetracene crystal, where there is a second exciton with 
a smaller oscillator strength, we observe a contrast de- 
pending on the orientation of the crystal of a factor of 
about 40. 
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